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ABSTRACT 

We have developed a universal approach to compute accurately the 
brightness of eclipsing binary systems during the transit of a planet in 
front of the stellar disk. This approach is uniform for all values of the 
system parameters and applicable to most limb-darkening laws used in 
astrophysics. In the cases of linear and quadratic limb-darkening laws 
we obtained analytical expressions for the light curve and its derivatives 
in terms of elementary functions, elliptic integrals and piecewise-defined 
function of one variable. In the cases of logarithmic and square root laws 
of limb darkening the flux and its derivatives were expressed in terms 
of integrals which can be efficiently computed using Gaussian quadrature 
formula, taking into account singularities of the integrand. 

Key words: stars, binary systems, eclipse. 



1 INTRODUCTION 

Recently several authors have developed 
algorith ms for the calculation of t r ansit l iRht curve s, 
see, e.g.. lMandelfcAgo]| (|2002l ).[Pail (|2008l ).|PaJI (|2012l ). 
However, the problem of calculation of the light curves 
is still relevant, because the existing algorithms are 
not applicable to all values of the system parameters 
for some limb-darkening laws. Besides, they do not 
allow sufficiently accurate calculations of the light 
curves for some limb-darkening laws. In addition, 
calculations of derivatives of the light-curve as a 
function of system parameters is important, because 
they can be used to solve the inverse problem of 
interpretation of the light curve . 

The paper iMandelfcAgoll (|2002l ) contains an 
analytical expression of the light curve by elliptic 
integrals, for the cases of linear and quadratic limb- 
darkening laws. In doing so, 13 variants of relations 
between the parameters are considered. For other 
limb-darkening laws (law of square-root and its 
power)only an approximate method of light-curve 
calculation at the radius of the planet more than 10 
times smaller than the radius of the star is being used. 
In this case, the accuracy is 2% of the depth of the 
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eclipse. In the paper |Paj (1201^)), directly, there is 
only an expression of the light curve in the linear and 
quadratic limb-darkening law, and the derivatives of 
the fight curve are calculated by difference methods 
(This work contains no direct analytical expressions 
for the derivatives), that is less favorable in terms of 
time and accuracy of the computation. In addition, 
none of the above works the logarithmic law of 
darkeningbare not considere d, which is the most 
prefer red f or early-type s tars (JKlinglesmithfcSobieskil 
l|l97d ) and lVanHamm3 (|l993D '). 

The approach presented in this paper allows us 
to calculate a light curve and with almost machine 
accuracy for any values of the parameters (including 
near singularities). Binary system parameters are the 
radii of the components, and the distance between 
the centers of the components in the projection 
on the picture plane. In general, the algorithm is 
uniform for all values of the system parameters, 
which significantly facilitates its implementation. We 
obtained analytical expressions for the transit light 
curve of the eclipsing binary system and for its 
derivatives in the cases of the linear and quadratic 
limb-darkening laws. These quantities are expressed in 
terms of piecewise-defined function of one variable and 
incomplete elliptic integrals, which can be compu ted 
with effective methods proposed bv lCarlsonI (ll99J). In 
the cases of the logarithmic limb-darkening law and 
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the square-root limb- darkening law the light function 
is expressed through integrals that can be efficiently 
computed using Gaussian quadrature formula. In this 
respect, the computation time of the light curve is not 
much more than the computation time by analytical 
expressions. 



2 MODEL DESCRIPTION 

We considered the model of the eclipse of a 
spherically symmetric star with thin atmosphere 
by another spherical opaque component (the other 
spherical star or a spherical planet). 

Fig. [l] shows the geometry of the stellar disks in 
eclipse. 

The brightness at the point of the disk of the 
eclipsed star with polar coordinate p is given by: 



J{p) = J(0)/ ( f- 



Here J(0) is the brightness at the center of this stellar 
disk, 



/(r) = (l-/(M(r))), 



(1) 



/i(r) = Vl — r^ 

/(M) = ^Afe/fe(M), 

where functions fk are such that /fc(l) = 0, are defined 
by the law of limb-darkening in question, and At are 
the coefficients of limb-darkening. 

In this paper we consider the following frequently 
used limb-darkening laws: 

Linear limb- darkening law, for which /(/i) = 
= \lfi{,i) = Ai(l - ^i); 

Square law of limb-darkening, which is 
characterized by the presence of the term 
Aqfq{n) = Ag(l — /x)^ in the expression for /; 

Logarithmic limb-darkening law, which is 
characterized by the presence of the term AL/L(/i) ~ 
— — Ai/iln/i in the expression for /; 

Square root limb-darkening law, which is 
characterized by the presence of the term AQ/Q(/i) — 
= A(3(l — -^//I) in the expression for /. TaKJKe nojiyuen- 
Hbie flajiee pesyjitTaTBi mojkho OMeBUflHBiM o6pa30M 
o5o5n],HTB na cjiynaii saKOua uoTeMHeHua k Kparo, xa- 



paKTepnsyeMOro ujichom v^ b BbipajKenun fljia ap- 
KOCTH, rfle I - nojiojKHTejiBHoe u,ejioe uhcjio. 



3 GENERAL INTEGRAL FORMULA FOR 
THE FLUX 

The decrease of the flux of the binary system due 
to eclipse is: 



L'^-L{D,R,,R^) = AL{D,R,,R^)^ // J(!R!)dR, 

S{D) 

(2) 
where L^ is the unobscured flux of the binary system, 
L is the obscured flux of the binary system, i.e. the 
light-curve value, S{D) is the area of overlapping 
disks, R is radius-vector of the point on the stellar 
disk. 

To calculate the integral ([2]) we introduce the 
functions: 



iTT, X < —1 

arccosa::, —I ^ a: ^ 1 
0, X > 1 . 



(3) 



and 



Then 



QX: 



dx 



V^, X ^ 
0, X <0 . 



= Q 



(4) 



(5) 



The relation ([5| is obtained naturally by noting 
that Az = Rearccosz, Qx = Rey^z for complex 
number z with Imz = and for the functions of 
complex argument arccos and ^ that are analytic 
continuations of the inverse cosine and square root 
of a real argument. Analyticity region is such that 
— TT < arg z ^ TT for each z. 

In the polar coordinate system the region of 
integration S{D) is given by: 



SiD) = <^ 



( p<R, 

n < ip i^ n 



p^ + D^ - Rt 
2pD 



(6) 



^ cos Ifi . 



In case of integration ([2| with respect to 
coordinate ip, for the values of p, which satisfy 



p'^ + D^ - rI 



2pD 
that 

p^ -HD^ ~ rI 

2pD 

Hence, the 



^ I variable ip takes the values such 



^ cos ip <^ \tp\ ^ arccos 



integration 



from 



p' + D' 



RI 



2pD 



p' + D- Rt 



2pD 

For the values of p, for which 
inequality 



p^ + D^ 



p' + D^^ Rl 
2pD 

tp is 

to 



r: 



2pD 



p^ + D' 



r: 



< COS ip 



2pD 
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Figure 1. The model of eclipsing binary system. The projection on the picture plane. Here the smaller component is a star 
or an exoplanet. The geometry of stellar disks in eclipse. Here ij, is the radius of the eclipsed star, Ro is the radius of the 
eclipsing component, D is the distance between the centers of the disks of the components, p, ^ are respectively the polar 
radius and the polar angle of a point on the disk of the eclipsed star. The origin is located at the center of the eclipsed 
star, the polar angle is measured counterclockwise from the radius- vector connecting centers of the star and the transiting 
component. 



holds for every value of if. In this case, integration 
with respect to ip runs from — tt to n. 



For the values of p, for which 



p"^ + D^ - rI 
2pD 



>1, 



the last inequality ((6]) is not satisfied for any values 
of ifi. Formally, at these values of p both integration 
limits by ip are set equal to zero. 

Next, using the notation Q and introducing the 
function 



integral in ^ can be rewritten as: 

R, <i'{D,p,R^) 

AL{D,R,,Ro)^ I pdp I dipj{p) = 

= 2J(0) I p'i'iD, p, i?„) / (^-|-) dp = 



= J{0)RlAL (^-|-, 1, 1^) = J{0)RlALiS,r) , (7) 
where r = -^ , 5 = -^ , and 



AL(5,r) = 2 / p*(5,p,r)/(p)dp = 



= / *(5,,/^,r)/(V^)rfp. (8) 



'i'{D,x,y)=A 



x^ + D' - y' 
2xD 



Note that AL(5, r) is the value of the decrease of the 
flux of the binary system when radius and brightness 



© 2012 RAS, MNRAS 000, [T]-?? 



M. K. Abubekerov and N. Yu. Gostev 



at the center of eclipsed star equals unity, radius of 
the second (eclipsing) component equals r and the 
distance between centers of disks equals 5. In view 
of Ul we can express the decrease of the flux as linear 
combination with limb-darkening coefficients: 

AL(<5,r) = ALo(<5,r) + A^AL,(^,r) + 
+A,AL,(^,r)+AtALL(5,r)+AQALQ(<5,r). (9) 

Unobscured flux L^ of the star is 
when R is the radius of the star, and 



I{^p)dp = 



= L({+A,l/+A,L^+AlL{+AqL^. (10) 

When both components of the binary system are the 
stars, unobscured flux L^ of the binary system is the 
sum of L^ for each star. For binary star and planet 
L^ equals L' for star. 

Let (7 be a function that g{p) is a separate linear 
term in the expression for I{y/p) (given by JU), 



g^ ' is one of its primitives: g{p) — 



dg^-^p) 
dp 



We 



consider the integral of the general form, which is a 
contribution to AL(5, r) caused by the term g{p) in 
the expression for I{^/p): 



AL,(5,r)=y ■^{5,^,r)gip)dp. (11) 



We note that for 5 > 0, r > 

lim ^((5, y/p, r) = TTQ(r - 5) , 

p—^O 

where 

fl, t>0 

[o, i<0. 
Using integration by parts, we obtain 

M.g{5, r) = *(<5, 1, r)g'--^\l) - 7re(r - %'"%)- 
\5^-r'-p)g^-'\p) 



X Q 



dp, (12) 



{p-{5-rr){iS + rf-p) 

where ((Sj is used for differentiating "if. 

For non-negative r and 5 the integrand in (|12p is 
non-zero only if (5 — r)^ < p < (5 -f r)^. Therefore 
the integral in (|12p is zero if \5 — r\ > 1. And if 
|(5 — r-| < 1, integrating can be performed over the 
interval ((5 — r)^,min ((5 + r)^, l)). In this interval 



S'^ + r^-p 
25r 



is a monotone function of p, so 



we can perform change of variable in integration in 
the following way: 



X — arccos 



r" — p 



25r 



Then 



p — S + r — 2rS cos x . 



(13) 



(14) 



Integration with respect to x will be perfomed 
over the interval 



0, arccos 



5^ + r^ -min((5-|-r)^,l) 
2l7 



Taking into account the fact that 



d fS^ + r^-p 

Tp '"''''''[ 2Sr 



and 



^ip-iS-rf)iiS + rf-p)^ 



S^ +r^ -min{{S + rf,l) 



= A 



2Sr 

5^ + r-2 - 1 



25) 



= *(5,r,l). 



we obtain: 



ALg(5,r) = *(5, l,r)g'"''(l) - 7re(r- - 5)g^-'-\Q) + 



*(5,r,l) 



/ 



(r — r5cosx))g^ ^'{5 +r — 2rJcosa;) 
(^2 _|_ j,2 _ 2r5 cos X 



dx 



(15) 



Expression (|15p is obtained assuming |5 — r| < 1. 
If |5 — rj ^ 1 value of ^(5, r, 1) = and integral in 
(|15|) vanishes. As noted above, when |5 — r| ^ 1, the 
integral in p2|l is zero because the integrand vanishes. 
Thus, the expression psp is valid for all positive values 
of (5 and r. 

By differentiating in pip integrand with respect 
to & and r we similarly find an expression for the 
corresponding partial derivative ALg: 



dS 



*(5,r,l) 



^-2./cos..(.^..^-2..cos.).. 



(16) 



dr 



*(5,r,l) 



2r g{5 +r — 2r6cosx)dx . (17) 



The contribution to the L-* caused by the term 
7(p) in the expression for I{yfp): 

L^=^(ff*~''(l)-5*"''(0)). (18) 
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4 INDIVIDUAL LAWS OF LIMB 
DARKENING. 

The expression for the decrease of the flux due 
to eclipse of the stellar disk with uniform brightness 
(with zero coefScients of limb darkening) can be 
obtained if we put in ^T^, (fT^ and {TTJ g{x) = 1, 
g~^{x) = X. Then: 



Lo = TT , 



ALo(5,r) = *(5,l,r) + *(5,r,l)r^--Q(5,r) , (19) 



aALo(5,r-) _ Q{5,r) 



and 



dS 



^ALo(<5,r) 
dr 



= 2*(5,r, l)r. 



(20) 



(21) 



Here 



^{S,T) = Q{{l-{S~Tf){{5 + rf-l)) 



Putting g[x) = ^J.{^/x) = ^/l-x, g' \x) = 
'2(1 — 1)2 , we get: 



L{ 



27r 



27r, 



ALi(5,r) = — e(r-5)+ 



+ Q 



l-{r-Sf 



2(5 + ; 



3(5 - r) 



n- 



■-{^5^-^) + {l-(r-5f){{r + 5f-l))F 



+ 



+ -^Q{l-{r-5f){7r^+5^-A)E. 



(22) 



Here 



n = n - 



4(5r V['(5, r, 1) 4(5r 



F = F 



E = E 



(r - 5)2 ' 2 1 - (r - 5)2 

*(5,r, 1) 45r 



2 1 - (r - 5)2 

1'(5,r, 1) 45r 

2 1 - (r - 5)2 



where F, E and H are incomplete elliptic integrals of 
the first, second and third kind: 



F{^\m) = 



de 



sjl-msin\e) ' 



E{d>\m) 






n{n;<l>\m)= f — 

J I 1 



'(0))v/l-msin2(e) 



The efficient algorithms for their calculations were 
suggested by ICarlsonI (j 19941 ). When \5 — r\ — 5> 
or |5 — r| — >■ 1, the limit of the term containing 
the factor fl in (|22|) is equal to zer o. Note th at a 
similar expression was obtained by IPall (120121 ) for 
the integral (primitive) of the appropriately chosen 
vector field along the limb of the eclipsed component. 
However, application of this expression for calculation 
of the flux of the system requires further account 
of its singularities. Expressions (|22|) and (I19|) give 
direct algorithm for calculating of the brightness, 
and the possible singularities are taken into account 
automatically by piecewise smooth functions of one 
variable ^ and Q 

aALi(5,r) _ 



dS 
*(5,r,l) 



'2r / cos(x)i/l — (5 — r)2 — 4 5r-sin2 ( - j 



dx 



35 



aALi(5,r) 
dr 
*(5,r-, 1 



Q{l-{r-5f) ((r + 5)^-l)F+ 



+ (l-5'' -r^)E 



(23) 



2r 



1 — (5 — r-)2 — 4 5r sin2 i-]dx 



= 'irQ{l-{r-Sf)E. (24) 
For term with linear limb-darkening coefficients in @: 

AL,(5,r) = ALi(5,r)-ALo(5,r), (25) 



t/ t/ 
^1 ~ ■L'O 



Assuming g{x) = x, g^ '{x) — x /2 we get 



T / TT 

^2 - T ' 



AL.(5,r) = -*(5,l,r) + 



+ y(25Vr^)*(5,r,l)- 

-i(5^ + 5r^ + l)Q(5,r) 



(26) 



The partial derivatives AL2: 
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aAL.(5,r)^ 25.2^(5,.,!) 



d5 



dr 



S'+r^ + l 
25 



Q{5,r) 



(27) 



2r (<5Vr^) 1'((5,r, l)-2rg(5,r), (28) 



Tf{k) 



-AiQk) 



cos 2x\'k — sin^ a- In (fc — sin x) dx 

(34) 



AL,((5,r-) = 2ALi(<5,r)-ALo((5,r)-AL2(5,r), (29) 



T / _ 97 / T / T -f - 



Further, we note that 

2 2 "^ V 4(5r 



Assummg (;(a;) — VI ^ xln(l — x), g^ {x) ~ 
= (1 - a;)^''^(4/9 - 2/31n(l - x)) we obtain for the 
logarithmic Umb-darkening law; 



L g . 



'-^±^^^H,5rf-^^i^+ 



dr 



+ 8rVr(5y^ ' ^ ^ 



4(5r 



(35) 



Vr[k) = 



2-^(Qfc) 



/ y^3-^ln(fc-sin^.) 



dx (36) 



Assuming (?(s) = ^1 ~ ^^ 9 {^) ~ iX~^) t^^ 

obtain the following expression for the case of square 
root limb-darkening law: 



ALi((5,r-) = ALi((5,r) {\n{A5r) - | 

2lT 

- — ln(45r)e(r-5)- 



8 /r 

-V5r 



,2yM('5-0^ l-(<5-r-) 



5r-y^ 



A&r A&r 

{5-rf l-{5-rY 
45r ' 4:5r 



(30) 



where 



47r 



ALi{S,r) = ^e(r-S)- 



8 4/. c 

- V4dr 
5 



2 0,0 



- (5ry! 



(5-r)2 l-(5-r)^ 



45r 
2 * ASr 



ASr 

2 1 /r „\2\ -\ 



QIJS- tY l-(5-r) 



45r 



, (37) 



^^(n.fc) 



-^(Qfc) 



/ 



(k — sin'^ x) ^ In (fc — sin'^ a;) 
n + sin'^ x 



dx, (31) 



where 



f--4(Qfc) 



T?(n,fc) = 






k — sin^ x) ■' , , 

^1— dx, (38) 



+ sm X 



7k{n,k) = 

l-A[Qk) ^ 

(k — sin^ x) ^ In (fc — sin'^ x) 
cos 2a:: -^ = dx 



n + sm X 



(32) 



The partial derivatives of AL^: 



i^-A{Qk) 



\n,k)= I 



cos 2a; 



/' 2 \ — 

[k — sin X) ■* 
n + sin^ X 



dx . (39) 



The partial derivatives of AL(; 



dAL3{S,r) 
dS 



- -^r<^^rS^f ('-^ 



(40) 



dALU5,r) ^^^^^^^^dAL,{S,r) 



as 



dS 



ASr 



(33) 



where 



Tf 



^^AiQk) 



(fc) = / cos 2a; (fc^ - sin^ a;) ^ dx (41) 
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and 



dAL3{S,r) 



dr 



4:r\'4:rS9r 



Qll^iS-rf 



4:ST 



(42) 



where 



'AiQk) 



9?{k) 



{k-smx)-'dx. (43) 



For term with square-root Hmb-darkening coefficients 
in ©: 



ALQ{S,r) = ALs{S,r) - AL„{5, 



(44) 



The formulas obtained for the square root, it 
is easy to generahze to the case hmb-darkening law 
contained in the expression for the brightness the term 
y7?, where I is an odd positive number, putting in 
(Hgi, (O add ini) g{x) = ^(1-x)', g^-''\x) = 



4 + « 



(1 



sl+l/4 



For an even I of non-multiple 



4, light curve and its derivative can be expressed by 
elliptical integrals similarly to as the formulas (|22p - 
(|24p were obtained. If / is divisible by 4 light curve 
and its derivative can be expressed by elementary 
functions similarly to as (|26p - (|28|) were obtained for 
quadratic limb darkening. 



logarithmic or fractional power singularity at f = 0. 
In the case of calculating of 9? vf , 99 , we put 7 = 1. 
For the limb darkening of the general form, which is 
characterized by the presence of the term \/fjJ in the 
expression for brightness (odd /) it is enough to put 
7 — I mod 4. 

By applying Gaussian quadrature formula, we 
can find the numerical value of the integrals with 
high precision, producing a relatively small number of 
elementary computations (the amount of computation 
of the integrand is proportional to required number of 
significant digits). But at the same time, an integrable 
function must satisfy certain conditions. In particular, 
this can be achieved if the higher derivatives of the 
integrand (or its non-singular component) is uniformly 
bounded on the section of integration. To reduce 
the computation of the integrals (|45p and (|46p to 
computation of the integrals that satisfy the above 

TV 

conditions, we divide the interval of integration (0, 

A{Qk)) sequence Xq > Xi > .... > Xm, such that 
Xo ^ ^ -- AiQk)), Xm = and 



k — sin^ XiJ^i 
k — sin^ Xi 



Xi — Xi+i ^ 



^ 2 for i ^ 1 and fc 7^ 1 , 



(47) 



max{r, 2} 



for all i < M andk. (48) 



In the case of (1 45 II we also require that the 
following inequality: 



n + cos Xi 
n + cosXi+i 



^ 2 for all i < M and k . 



(49) 



5 NUMERICAL CALCULATION OF 
INTEGRALS 

Thus the calculation of the brightness for the 
logarithmic and square-root limb-darkening law is 
reduced to the calculation of integrals J>f , y^, 9^ "Pg , 
Tr, Vi, V2, V? 7f, 7? (depending on parameters). 
These integrals can be represented in general form: 

|--4(Qfc) 

-. . f V(k,x)K{k -sin^x) , , 

7{n, k) = / ^ ' V^ ^ dx , (45) 

J n + sm X 



for yf, 92, ??, 72 or 



-AiQk) 



y(fc) = 



for y^ 7s, 7r, 7' 



V{k,x)K {k-sin^x) dx , (46) 
'^ Here V{k,x) = 






= "Y^ Ui{k)vi{x) where Vi{x) are some trigonometric 

polynomials, n > 0, fc > 0. We denote the maximum 
degree of these trigonometric polynomials as r. 
K{y) ~ ^/y^'^y or K(y) — ^pip , respectively, has a 



If A: > 1, inequality from (|47|l also holds for z = 0. 
If fc < 1, then 

k — sin^ X\ 



(Xo-X.)sin(2Xo)^2^"^^^^^°/2- ^^^^ 

(|47p - (|50p can be used as recurrent relation, 
allowing us to construct the sequence Xi. 
Thus, 



7{n,k)= ^ 7^{n,k), 



7{k) = Y. 5'»(^) ■ 



where 



Jfi 



- , , , V{k,x)K(k-w^^x) , , 

7^{n,k)= I !,,:., 2 ^ -dx, (51) 



n + sm x 



and 



^i+l 



Xi 



7i{k)= I V{k,x)K{k-s.in x) dx, (52) 

For fixed values of n and k the two last integrals can 
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be represented in general form: 

where U{x) = ^''^'7 for ^ and U{x) = V{k, x) 

n + sin X 

for ([5^. 

By linear substitution of variable of integration 
x{t) = Xi+i + t{X, - X,+i) 
in (|53[) . we turn to the integration from zero to unity: 



1 
9. = {X, - X,+i) I U{x{t))K [k - sin'^(a:(i))) dt . 



(53) 
In this form, J'i can be computed by applying the 
Gaussian quadrature formula: 

1. N 

h{t)uj{t)dt ^ ^ wih{ti) . (54) 

'=1 

Here u}{t) > OVt £ (0, 1), nodes f; are the roots of 
the polinomial HNit), where {Hi\ is the system of 
orhtogonal polynomials with weight uj in the interval 
(0,1): 



Hi{t)Hj{t)oj{t)dt = Ofor / / j . 



wi can be found as the solution of the system of A'' 
linear algebraic equations, which can be obtained if 
we put h{t) = 1, h(t) = i, . . . , h{t) = t^ in ^ and 
replace the approximate equality with exact equality. 
A'^ can be adjusted so as to ensure the required 
accuracy of calculation of Toiin, k),'J>ai{k) and can be 
the same for all values of i, n, k. N is of the same order 
of magnitude as the number of significant digits in the 
result, and this allows to calculate the integral with 
the required accuracy in a reasonable time. So, after 
calculation of the roots of polynomials xi and weights 
Wi (this may take a while), we can re-use them for 
computing 'J'oi{n,k),'J>ai{k) for all i,n,k. 

In the case of i > or of fc > 1 we 
put in (154|: uj{t) = IVi G (0,1), h{t) = 
= {X, - X,+i)U{x{t))K (fc - sin^'c(t)). Then 



7^ 



'^■wih{ti 



Note that in this case Hi{t) = Pi{2t — 1) where Pt are 
Legendre polynomials. 

In the case of i > 0, fc < 1 and logarithmic 
limb-darkening law {K{y) = ,/ylny) we represent the 
integrand from (|53p in the form: 



U{x{t))VT^ 



k-sm^{x{t)) 



X 



In 



1-i 

fc- sin^(a::(i)) 



4- ln(l - i) 



Next, we put in ((54|: Lu{t) = ^1 - tVt G (0, 1), 

'fc-sin^(a::(i)) 



, , -, r-r, , s^ Ik — sin^(2;(f)) , 
h{t) = U(x{t))\l -_ ^—^-^ In 



1-i 



Note that here Hi{t) = p/^'°'(2i - 1) where p/ = '°' 

N 

are Jacobi polynomials. Let Si = "^ wih{ti). 

1=1 
Next, we put in (fSi)) : uj(t) = 

= -y/T^^\n{l ~ tyit G (0, 1), 



hit) = -U{x{t)) 



fc-sin^(x-(i)) 

T~~t ■ 



The polynomials corresponding to this value of oj 
can be obtained through the standard procedure of 

JV 

orthogonalization. Let 82 — '^ wih{ti). Then Jo ~ 

iSi+S2){Xo-Xi). 

In the case of i > 0, fc = 1 and logarithmic 
limb-darkening law {K{y) — y^lny) we represent the 
integrand from (|53p in the form: 

ln(^2£M))Uln(l_t) 



2U{x{t))cos(x(t)) 



Next, we put in ((54|: aj(t) = IVi G (0, 1), 

h{t) = 2Uix(t)) cos(x(f)) In fSS^^M. 

N 

Let Si = Y:,wih(ti). 
1=1 
Next, we 



put 

= -ln(l-f)Vf G (0,1), 



uj{t) 



h(t) = ~2U{x{t))cos{x{t). 



Let 52 = E wih{ti). Then Jo ~ (5i + 52)(Xo - Xi). 

In the case of i > 0, fc < 1 and square-root limb- 
darkening law {K{y) = ^y^), we put in H54|) : u{t) — 
= (l_i)iVtG(0,l), 



h{t) = U{x{t)) 



fc — sin^ x{t) 



Note that here Hi{t) = P- * ■"' (2i- 1) where P^ ^ '"' are 

N 

Jacobi polynomials. Then CPo ~ {Xo — Xi) E wih{ti). 

In the case of z > 0, fc = 1 and square-root limb- 
darkening law, we put in (|54ll : Lj(t) — ^/l~^\/t G 

e(o,i). 



h{t) = U{x(t)) 



<t)) 



1-t 



Then'yo~{Xa^Xi)Y.wih{ti). 
1=1 
Calculations show that in all cases of the 
applications of the Gauss quadrature accuracy of 
19 significant decimal digits (corresponding to 80-bit 
machine numbers) can be achieved by choosing the 
N to be 14. Value sets of points ti and weights Wi 



© 2012 RAS, MNRAS 000, [T]-?? 



A universal approach to the calculation of the transit light curves 9 



corresponding to each of the considered forms of the 
function uj, can be downloaded from the Internet, 
along with other materials (see Conclusion). 



6 CONCLUSION 

We have derived the expression for the calculation 
of the eclipsing binary flux and its derivatives. 
We considered the linear limb-darkening law, the 
quadratic limb-darkening law, the logarithmic limb- 
darkening law, and the square root limb-darkening 
law. In general, the decrease of the flux is given 
by the expression @. In p9p - (|2ip . ALo corresponds 
to uniform brightness and its derivatives, it is 
expressed in terms of easily computed piecewise- 
defined functions of one variable ^ Q and Q Q. 
AL; corresponds to the linear limb-darkening law, 
given as linear combination of ALo and ALi (|44|) . 
where ALi with its derivatives is expressed in terms 
of incomplete elliptic integrals in (|22|l -(|24[). AL, 
corresponds to the quadratic limb-darkening law, 
given as linear combination of ALo, ALi and ALi 
(|29|) . where AL2 with its derivatives is given in 
(|26|) - (|28| ). ALi corresponds to the logarithmic limb- 
darkening and ALq corresponds to the square-root 
limb-darkening expressed by two- and one-parametric 
integrals. Further, we described how these integrals 
can be found numerically by multiple application of 
the Gaussian quadrature formula. It is important 
that the nodes for this formula can be found once 
and re-used for the calculations for different values 
of parameters. Also, the general integral form (|15p - 
(|17|l of the flux component allows us to extend this 
approach to other limb-darkening laws. 

The algorithm desc ri bed ab ove was tested by 
lAbubekerov et all (|201(]| . I2OIII '): iGostevI (|201l| ) for 
the interpretation of the high-precision polychrome 
light-curves o f the binary s ystem with exoplanets 
HD 209458 iBrown et all I2001I ) , HD 189733 



IPont et al.l (|2007|) and monochr ome light-curves o f 
Kepler-5b, Kepler-6b , Kepler- 7b |Koch et all f{20ldi ): 
iDunham et all (|201(]| '): lLatham et all (|201(]| '). 

The algorithm is implemented in ANSI C 
in the form of the functions for computation 
of the individual component AL((5, r) and its 
derivatives. This implementation is available from 
http://lnfml.sai.msu.su/~ngostev/algorithm.html 

This work was supported by the President 
of Russian Federation grant MK-893.2012.2, RFBR 
grant 12-02-31466. 
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